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£> . Abstract. A new computational method for solving the nucleon-deuteron breakup 

scattering problem has been applied to study the inelastic neutron-deuteron scattering 
on the basis of the configuration-space Faddeev equations. This method is based on the 
spline-decomposition in the angular variable and on a generalization of the Numerov 
_u. ' method for the hyperradius. The Merkuriev-Gignoux-Laverne approach has been 

f^ . generalized for arbitrary nucleon-nucleon potentials and with an arbitrary number of 

partial waves. Neutron-deuteron observables at the incident nucleon energy 14.1 MeV 
have been calculated using the charge-independent AV14 nucleon-nucleon potential. 
Results have been compared with those of other authors and with experimental 
S^ ' neutron-deuteron scattering data. 



1. Introduction 

The last two decades brought tremendous progress in both theoretical and experimental 
study of nucleon-deuteron scattering. More accurate experimental data were compared 
with rigorous calculations in the framework of the Faddeev equations with high- 
precision nucleon-nucleon potentials also including model three-nucleon forces [U [2]. 
Nucleon-nucleon (NN) potentials used in rigorous three-nucleon (3N) calculations are 
charge-independent AV14 [3], AV18 [I], CD-Bonn [5] and several Nijmegen potentials 
[6]. Among three-nucleon forces (3NF) are Tucson-Melbourne (TM) and its various 
modifications [7], and Urbana potentials [H]- Based on the chiral effective field theory 
(EFT) NN and 3N potentials have been constructed [9] and used in rigorous 3N 
calculations [2]. 

In spite of these enormous achievements in the 3N studies, there are several 
important cases where the rigorous 3N calculations have failed to explain the data 
[TO] and these discrepancies are established with very high precision. Among the most 
important discrepancies are the A y puzzle in nucleon-deuteron (Nd) elastic scattering 
|llj . the star configuration in the Nd breakup reaction [12] and the quasi- free scattering 
(QFS) cross section pp. 

Recent calculations of the proton-deuteron breakup cross-section, and tensor 
analyzing power data for a symmetric constant relative energy (SCRE) geometry 
reported in [13] revealed a serious disagreement between the theory and experimental 
data. Inclusion of the TM or Urbana 3NF did not solve the problem. The chiral 
EFT showed the same deviations from the data. Including the Coulomb force partially 
improved the situation for the cross section. On the other hand predictions in [H] with 
the 3NF forces and Coulomb interaction taken into account showed a good agreement 
with the data. At the same time, the recent study of Nd breakup in symmetric 
star configurations based on Faddeev calculations with CD-Bonn NN potential [15] 
showed significant discrepancy for nn QFS configuration and a good agreement with 
cross-section data in the case of pn QFS breakup geometry. Thus the situation with 
theoretical predictions for the 3N observables remains still far from completely resolved. 
This motivates us to continue searching for new computational procedures and input 
potentials, which may allow to overcome the mentioned disagreements. 

In this article we present an alternative method for the study of the neutron- 
deuteron (nd) system based on the direct numerical solution of the Faddeev equations 
in configuration space. This approach was initiated by Merkuriev et al. (MGL) [16] 
who derived general formulae for nd breakup scattering. This method has been initially 
applied to the study of nd and pd elastic and breakup scattering with nuclear interaction 
limited only to S-waves and with simple NN potentials [17]. In the present work we 
generalize the MGL approach to any high precision realistic potentials for nd inelastic 
processes. 

This paper is organized as follows: in section 2 we describe calculations in 
configuration space starting with the general formalism in subsection 2.1 followed by 



Numerov method in subsection 2.2. Our novel method for solution is given in subsection 
2.3. Formulae for calculating elastic and breakup amplitudes in MGL basis are presented 
in section 2.4. Comparison of our results with the previous calculations and with the 
data are discussed in section 3. Finally, discussion and conclusion are given in sections 
4 and 5. 

2. Three-nucleon Faddeev calculation in configuration space - our new 
computational method 

2.1. Formalism 

The starting point for studying interactions between nucleons in three-body systems is 
the solution of the Schrodinger equation H \1/ = E^fr for nuclear Hamiltonian 

tf = -^Ev? + E^ (+ E *$«)> (i) 

z '" L j=l j<k j<k<l 

where Vjk and Vjki are the two- and three-nuclear potentials, respectively. In this study 
we neglected three-nucleon forces Vju- 
Writing the total wave function as 

^ = $ 1 + $ 2 + $3 = (l + P+ + p-)$ li (2) 

the Schrodinger equation for three identical particles can be reduced into a single 
Faddeev equation, which in Jacobi's vectors X\, y\ has the form 

(A* 1 +AftW(£i)-.El$(£i,jfi) = -V(x 1 )(P + +p-)$(xi,yi),(3) 

where the operators P ± are the cyclic permutation operators for the three particles which 
interchange any pair of nucleons (P + : 123 — » 231, P~ : 123 — > 312), and Ti 2 /m=AlA7 
MeV-fm 2 . As independent coordinates, we take the Jacobi vectors x a) y a . For the pair 
a=l, they are related to particle coordinates by the formulas: 

xi = r 2 -r 3 , yi = — n, (4) 

for a=2,3 one has to make cyclic permutations of the indexes in Eq.(5). The Jacobi 
vectors with different a's are linearly related by the orthogonal transformation 

{i)={-Z %){%)' ca - + *= i - (s) 

where 



To perform numerical calculations for arbitrary nuclear potential, we use MGL 
approach [IB]. To make the angular analysis of Eq.(j3]), we use a basis proposed in [16J. 
This basis is intermediate between LS and Jj coupling schemes. 
Let: a, 1 and J = a + 1 be the spin, orbital and total angular momenta of the pair 



23. s = 1/2 + J the total "spin" of the system 123 considering the pair 23 as a particle 
with "spin" J. A the orbital momentum conjugate to y, that is of the relative motion of 
particle 1 respective to the cm of pair 23. M = A + s the total angular momentum with 
its ^-projection M z . To this we have to add the isospin part. If the isospin of the pair 
23 is t, then the total isospin is T = 1/2 + t with its projection T z . Since in the nd case 
T = 1/2 it need not be shown explicitly. The set of quantum numbers {As/a J} = a 
defines a state in this basis. 

Correspondingly in this basis the spin-angular-isospin eigenf unctions have the form 



Z a {% y) =< x, y|a >= W\y) <g> [ X l/2 ® 3^(x)] 



MMz 



T,T Z 
%/2,t 



(7) 



For nd scattering the Faddeev equations for partial components can be written in 
the following form (here we omit the index 1): 



E + g(«9 x 2 + d 2 y ) - v%(x, y)] $*» so > M o(x, y) 






iAo,so,Mo / 



=,Ao,so,Mo 1 



^°' so ' m °(x : y) + du^2g fiy (y/x,u)^ '"{x,y 

J-i 7 



(8) 



The geometrical function g^(x,y,u) is the representative of the permutation 
operator P + + P~ in MGL basis [16J: 

g a ' a {y/x,u) = g a > a {0,u) = g a > a (9,9') 



(-1) 



X+X'+J+J' 



[(2 J + 1)(2 J' + l)(2s + l)(2s' + I)} 1 ' 2 J2(2S + 1)(2L + 1) 

LS 



(9) 



j I a J 1 f V a' J' 1 f A I L )( X' l' L ) 
\ 1/2 s S J [ 1/2 s' 5 J j 5 M s J j 5 M s' J 

< Xl/2a'Vl/2,t'\ P+ \xf/2aVl/2,t > h\>l'Xl(y/x, «) ■ 

Function ft, is the representative of the permutation operator P + + P~ in the 
A + I = L basis: 

g7(-iy +L (2A+ 9 1 A K f + 1) [(2A)!(2/)!(2A- + 1)(2/' + 1)]^ 



x'y' ' 2 A 

£(-l) fc (2A; + l)P fc 

(Vs) X2+h 



u 



k=0 



„,Ai+Zi™A2+i2 

x 1 +x 2 =x,i 1 +i 2 =i y x 



[(2A 1 )!(2A 2 )!(2/ 1 )!(2/ 2 )!] 1 / 



_£(2A" + 1)(2Z" + 1) 
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(10) 



The index k runs from zero to (A' + V + A + l)/2. The (...) are the 3j symbols: 

1 



/ 3\ 32 33 
\ mi m2 m3 
The centrifugal potential is 
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(ii) 



and nucleon-nucleon potentials are v aa >(x) =< a\v(x)\a' >= 5\\i S ss i5 aa >6jj'Vf/ , where 
vf/ are the potential representatives in the two-body basis y i(T z (x) (most often 
abbreviated as 2cr+1 lj). 

The set of partial differential equation Eqs.([HD m ust be solved for functions 
satisfying the regularity conditions 

$M /o (0, 0) = $W/ 0( ^ Q ) = QXosoMofa tt/2) = (12) 

The asymptotic conditions for nd breakup scattering has the following form |18j : 

qXosoMo^ y) „ | [5 XXo 5 sso 5 al 5 n Uqy) + ( - Hqy) + ua(^))4X J<M^) 



(13) 



X 2 = x 2 + y\ K 2 



mE 

— 2~, x finite, y — > oo, 



where tpi is I — th component of deuteron wave function (/ =0,2), and j and y are the 
regularized spherical Bessel functions. 

The matrix of partial elastic amplitudes a has the structure 
r]exp(2i5) — 1 



a XsX o s - 2i ' ^ ' 

where rj and 5 are the inelasticities and scattering phases. In Eq. (ll3j) the amplitudes 

A^ x are related with physical breakup amplitudes. The inelastic scattering 
amplitudes are given by the sum 

<U(*) = <W') + f d^9a^,n)AfXj0')- (15) 

To simplify the numerical solution of Eqs.([E|), we write down Eqs.flE)) in the polar 
coordinate system (p 2 = x 2 + y 2 and tan# = y/x): 



^ + ^^ + ^ + ^-^^ f/ " S0M °^^ 



(16) 



Y / v af} (p,9)[u^ Mo (p,9)+ / duJ2g^(e,u,e'(e,u))u^ soMo ( P ,e') 

Here the first derivative in the radius is eliminated by the substitution U = p _1 ' 2 $. In 
Eqs. (15-16) the angular variable 9' is defined by 



1 P\ ^ 

cos 9 (u, 9) = — cos 9 cos 9 sin 9 ■ u + - sin 9. 

y ' 4 2 4 



;i7) 



2.2. Numerov method 

Modification of the Numerov method for the set of the differential equations ( TT6l) does 
not present any difficulties in principle. As is well known, the Numerov method is an 
efficient algorithm for solving second-order differential equations. The important feature 
of the equations for the application of Numerov's method is that the first derivative has 
to be absent. The aim of this method is to improve the accuracy of the finite-difference 
approximation for the second derivative. Starting from the Taylor expansion truncated 
after the sixth derivative for two points adjacent to x n , that is for x n _i and x n+ \ one 
sums these two expansions to give a new computational formula that includes the fourth 
derivative. This derivative can be found by straightforward differentiation of the second 
derivative from the initial second-order differential equation (see the details in [T9]). 
For brevity, we omit the corresponding derivation and present only the final formula of 
Numerov's method for Eqs. (TT6l) . not indicating indices A s M : 

r 12 2ApT a (9h 

-F+(A^ + (1 + — )-^"J^(Pi-i^) = 

$><* (Pi. 9) - Apv' a0 ( Pj , eMUpfa-!, 6) + Y: [ + _ de'g M {6, ^ 7 (p,-i, &)) 



+J2(iov a ^ Pj ,e) + (Ap) 2 v^(p j ,e))(Up(p j ,e) + J2j e _ de'g^o^u^e')) 

-[E + j^ + (i - —)-^\U a (p j+1 ,0) 

+ Y,(v a ^p J ,e) + Apv' aP (p J ,e))(u^ +u e) + Y,rde'g^(e,e')u 1 (p J+1 ,e')) = o, 

a n, J 6 



where 

T (R\ - ?L - ^ + 1 ) _ A ( A + 1 ) , I 
a{) dO 2 cos 2 ^ sin 2 ^ + 4' 

In Eq. ( TTHl) pj is the j — th current point for hyperradius p G (0, R ma x) in the radial 

grid (j = 1, 2, . . . , N p ), Apj is the radial step-interval. 

To ensure the accuracy of order (A#) 4 for the approximation in the angular variable, 

Hermitian splines of the fifth degree have been used (see Ref. [20]). These splines are 

local and each spline S a i(x) is defined for x belonging to two adjacent sub intervals 

[xj_i,Xj] and [xj,Xj+i]. Their analytical form is fixed by the following smoothness 

conditions: 

S oi {x i - 1 ) = 0, S ai {x i+1 ) = 0, a = 0,1, 2, (19) 

and 

Soi(xi) = 1, S' 0i (xi) = 0, S$i(xi) = 0, 

Su(xi)=0, S' ll (x i ) = l, SUxi) = 0, (20) 

S 2 i(xi) = 0, S' 2i (xi) = 0, S^Xi) = 1. 



Expansion of the Faddeev component into basis of the Hermitian splines has the 
following form: 

2 Ng+1 

u a ( P ,9) = Y: E 5 ff #)c;.( P ), (21) 

a=0 j=0 

where N$ + 1 is the number of internal subintervals for the angular variable 9 G [0, vr/2]. 
To reduce the resulting equation (ITS]) to an algebraic problem, one should explicitly 
calculate the derivatives of NN potentials v a/ 3(p,9) with respect to p and the second 
derivates of splines S a j(9) with respect to 9. It is convenient to express the second 
derivative of component U a with respect to 9 through U a itself using Eq. (j2T|) . Upon 
substituting the spline expansion ( 1211) and expression for its second derivative into 
Eqs. (TT8|) . we use a collocation procedure with three Gaussian quadrature points per 
subinterval. As the number of internal breakpoints for angular variable 9 is equal to 
No, the basis of quintic splines consists of 3A^ + 6 functions. Three of them should be 
excluded using the last two regularity conditions from ( fl2l) and continuity of the first 
derivative in 9 of the Faddeev component at either 9 = or 9 = tt/2, as the collocation 
procedure yields 3N e + 3 equations. Finally Eqs. (TT8l) for the Faddeev components are 
to be written as the following matrix equation: 

A 1 U 1 + G 1 U 2 = 0, 
BjUj^ + A,U J +G J U J+1 =0, j = 2,...N p -l, (22) 

Bn p Un p -i + A Np U Np = -G Np U Np+1 . 

In this equation vector Uk = U{pk) has dimension N in and matrices A, B, G, derived 
from Eq.(18), have dimension N n x N n where N n = N a x N c , and N a is the number 
of partial waves and N c = 3Ng + 3 is the number of collocation points in the angular 
variable 9. 

2.3. The novel method of solution 

To derive equations for calculation of breakup nd amplitudes, the method of partial 
inversion [19] has been applied. We write down Eq.( !22l) in a matrix form: 

(D*U) i = -6 iNp G Np U Nl)+1 . (23) 

Here matrix D is of dimension N p N in x N p N in , and N p is the number of breakpoints in 
the hyperradius p. The form of this equation results from keeping the incoming wave in 
the asymptotic conditions f|T3|) . As a consequence, the right hand part of Eq. (122]) has 
a single nonzero term marked with index N p + 1. Sparse (tri-block-diagonal) structure 
of matrix D optimizes considerably the inversion problem. Hyperradius p n+ i = R ma x, 
where R max is the cutoff radius at which the asymptotic conditions Eqs. (!T3|) are 
implemented. By formal inversion of the matrix D in Eq. (23), the solution of the 
problem may be written in the following form: 

U 3 = -Dj r lG Np U Np+1 , j = 1, 2....N p . (24) 



In Eq. ( |24p one should consider the last two components of vector U: 

Un p -i = -D Np _ lN G Np U Np+1 
U N = -DZ N G N U N+1 . 



(25) 



Provided R ma x is large enough, the vectors Un p -i, Un p on the left side of Eqs.(25) may 
be replaced by the corresponding vectors obtained by evaluating Eqs. (fTBl at the radii 
P — Pn p -i an d p = Pn • As a result we obtain a set of linear equations for the unknown 
amplitudes a and A: 

a ■ v Np „! + mxp-! ■ A = JF Np -i ^a) 

a ■ v Np + m Np ■ A = F Np . 

For the sake of brevity, we do not display here the explicit form of vectors Vj , J-j and 
matrices rrij. As R ma x — > oo the set of equations (1261) has a constant a as a solution. At 
finite R m ax its solution is a vector a with generally different components corresponding to 
different angles. We follow the method of S.P. Merkuriev [18J, which consists in selecting 
the components of a in the region of the maximum of the deuteron wave function, where 
a turns out to be independent of the angle. 

Furthermore, we propose a new method for a more adequate calculation of the 
amplitudes. The set of linear equations ( |26l) is over determined, since the number of 
equations is 2 • N in and the number of unknowns is iVj n +l. Therefore it is natural to use 
the least-squares method (LSM). One can apply it by two ways. In the first one, it is 
needed to express the breakup amplitude A from the lower equation (126J1 and substitute 
it into the upper one. As a result one has the following expression: 

o-v = F, (27) 

where vectors are defined as follows: v = v n -i — tn n ^ 1 m^ 1 v n , F = J-'n-i — rn n -i m' 1 Tn- 
According to LSM one should to minimize the following functional 

\\a- v-F|| 2 = min. (28) 

Differentiating this expression in Re a and Im a we obtain 
(v*,F) 



(29) 

where (£*, /) is an ordinary scalar product. 

In the second way, it is needed to express the elastic amplitude a from the lower 
equation (1261) using the scalar product: 

(v*,T n - m n A) 



v V n> V n) 



(30) 



Substituting a from Eq. ( 130]) into the equation ( )26|) . leads to the set of linear equations 
m n -iA-v n -i— r - — = Tn-! - u„_i— r. (31) 



The explicit form of Eq. ( 131 j) is as follows 

l2j=l \jn n _ 1 ^j — (v* <v ^ Hk=l V n,k m n,kj]Aj 



(32) 



T, 



n-U Un-l,i \ v * >Vn ) ) * 



1 TV 

J. , ..., iV c . 



Solving the set in Eq. (152]) . we get the breakup amplitude A Substituting the obtained 
breakup amplitude into Eq. ( 150]) . one may compute the elastic amplitude a. Note that 
one can apply Eq. (130]) to calculate the elastic amplitude a either in the components 
or via a scalar product. In the first case, the components of a are practically equal to 
a constant for all angles 9 G (0,7r/2) and this constant coincides with the value of a 
calculated by using the scalar product to the fourth decimal. It should also be noted 
that the elastic amplitudes calculated by the method from [18] and LSM coincide with 
this constant to the same accuracy. To control the accuracy of calculations, all methods 
are used. 



2-4- Observables 

To calculate observables for elastic scattering of nucleon from deuteron in the direction 
q' (initial direction q is along the z-axis), one has to derive the equation for the elastic 
amplitude as a function of scattering angle. Omitting this derivation, we represent the 
final expression for this amplitude in MGL basis: 



a<j'J'a z ,J z 



#) 



nnn 



•A-A' 



syMM z 
W'M Z - 



M X's' As 
is'<K.+J* rysa z +Jz 



f 2A + l 

An 



(33) 



-J's'a'+J. 



s~iMM z /~ii <"o" z -t-j " z ris(j z +j z y 

,Ly \0,sa z +.J z { -'l/2a' z ,lJ z ^l/2<T,lJ z I XM z 



-J z (q') a A's'As> 



with M z = a z + J z . 

In Eq. (133]) (x'cr^a, crj and J'J' Z (JJ Z ) are spin and its projection for incoming 
(scattered) nucleon, and the deuteron in the rest (scattered deuteron), respectively. 
Thus, the nuclear part of the elastic amplitude is a (2 x 2) £g> (3 x 3) matrix in the spin 
states of nucleon and deuteron, depending on the spherical angles 9 and 0. 

The spin elastic observable formulas can be taken from the review [I]. They are 
expressed via spin 2x2 matrices Oi for the nucleon and 3x3 matrices V% and V%k for 
the deuteron. The latter are related to the deuteron spin matrices Sf. 



S 



( 1 \ 

1 1 

V° i °/ 

Si, Vik = 3/2(SiSk 



1 

V2 



o n 



1 
V2 



One has Vi 

3(>J a ;D a ; OyOy). 

Nucleon analyzing powers A^ are 
Tr [aakO^) 



SbSi) 



( -i \ 

i — i , S z 
v i i 

I zz = <JQ Z Q Z Zl , 





(l 













\ 


-1 


an 


d V xx - 


- V 

1 '1 



\ 



.(34) 



yy 



A k 



Tr (aa^) 



(35) 
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If the scattering plane is the xy plane and the y axis points to the direction qxq' then 
due to parity conservation A x = A z = and the only non-zero component is A y . The 
deuteron vector and tensor analyzing powers are defined as 

At= iimp, Aj ^^mgi. (36) 

Tr(aat) ' Jk Tr(aat) v ; 

Parity conservation puts A X ,A Z , A xy and A yz to zero. So the non- vanishing and 
independent analyzing powers are defined by 

iTu = ~Y^v ^"20 = —teA-zz-, T 2 \ = — -t=A xz , T 22 = — -y=(A xx -A yy ).(37) 

Also spin transfer coefficients are given in the review [lj. They have the same 
structure as the quantities above, with slightly different matrices to be inserted between 
a and a*. 

In the case of Nd breakup scattering expression for physical breakup amplitude to 
calculate breakup observables has much more complex form. Below we present the main 
details of its derivation. The asymptotic of the wave function with the given incident 
plane wave is defined as follows 

X^T ^ Aoz ' SOz {q) ^ sin 9 cos 9 Z ^ y ' xj ' 

where the initial plane wave has the form 

cttf(q) = ^X^oJ^ AoJ*(q), (39) 

and Z Q (x, y) are the spin-angular-isospin eigenf unctions. 

The exponential factor is multiplied by the function depending only on angles, that 
are directions of vectors q, x and y . This gives a probability to find the breakup particles 
at given directions and so with the amplitude for the breakup with given directions of 
p' and q'. So the breakup amplitude is 

4-7T _ _ A^oso,Mo(Q\ 

A,o,Up, q') = - E ct!°;f°(q) E ? e Je Za ^ &> (40) 

y A ,A Oz ,Mo,M Oz a Sill (7 LOS (7 

where now 9 = arctan(g'/V) 

The only thing necessary to pass to the formula for the breakup is to project this 
amplitude onto the state with given projections of spins of the three particles fix, /i 2 , fi 3 . 
This will obviously be given by introducing into the sum over a in Eq. (HO"]) the projection 

<°Z% 3 (f>'^) =< AW 3 |^(q', p') > • (41) 

As the result we get the breakup scattering amplitude as a function of final nucleon 
momenta in the following form 

A(p, q, HlTi[l2T 2 fl 3 T 3 \q, s , s 0z ) 
— V^ jMqMoz (fJ A/N A*,Aq,sq \y ) XqsqMq / */n ( ac)\ 

- l. VnM 2 ^ ST -3,«ip>qj sin ^ cos ^ c Ao z , S o z iqj ^> 

^ a,7r,Ao,Aio,Aoz,JVJoz 
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where //j and T{ i = 1,2,3 are the spin and isospin projections of the three nucleons, 
9' = arctan(g'/p'), and summation goes over A 0z + s 0z = M 0z . A a ^ so is the spherical 
inelastic amplitude defined in Eq. ( Tl5l) . The ^-coefficients are 

CSWs^p'^O = (-i) A+J+M - 1/2 [(2J + i)(2, + i)] l/2 C7^ 3 

1 1 



Ci T2+ r 3 CI F , ^ V[(2L + l)(2S + l)] 1/2 C? Sz , Ctit° z s ss 

yjfo.S'^tf) (43) 



2 ' z >2 '<> 2 2 ' z '<»■".« ' ' a — J' 

/ a J ) I A I L \ ^ 7 l,m 0z -Sz/£.i ~/\ 



[1/2 s S j { S M s 

where S z = [i\ + /x 2 + /i 3 and t=0 or 1 according to antisymmetry condition I + a + t 
odd. 

The breakup differential cross-section is then 

dSj/tfq = 3^l^(P / '^^ r i^r2^3r3|q,So,Sa B )r (44) 

Note that d 3 p'q'd 2 q' = p'd 2 p'd 3 q' 

This formula may be transformed to the form which is used by the experimentalists 
in the lab. system. 

|^(p,q,/iiri/i2r2/i3r 3 |q,So,So 2 )| (45) 



d^ 2 kid 2 k 2 gif 3 V£ 
where 

D = A; 2 (2A; 2 - k 2 (k /ab - k x )) 2 + fc*(2A;i - k^k^ - k 2 ))' (46) 

and S is the arc lenght along the allowed curve in the E\ — E 2 plane: 



dS = dEi r 



2 

k 2 2kx -ki(k /a6 -k 2 )\ 



(47) 



"\| \k 1 2k 2 -k 2 (k lab -k 1 )J 

This cross-section refers to the breakup experiment in which all individual spins 
and isospins of the final nucleons and the spin of the initial neutron-deuteron system are 
given. The unpolarized cross-section with the isospin projections of the final nucleons 
given is obtained by summing Eq. ( 145]) over /ii,/i 2 ,yU 3 and averaging over s and its 
projections So^: 

d 5 a 
dSd 2 kid 2 ~k 2 (AQ \ 

2^ |-4(P,q,Atin/W 2 r 2 /i3r 3 |q,So,So 2 )| 



qK v J-J ° Ati,^2,M3>so,so z 
To write this expression in a more convenient form, one can introduce a matrix 
between initial and final spins: 

*4(p', q'lqWa^sKso* =^(p / ,q / ^in^2r2/X3r 3 |q,so,So z ) (49) 
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In terms of this matrix the unpolarised cross-section can be written as 



fa _ ^ mk ^l Tl 0A} (50) 



dSd 2 k t d 2 k 2 qK 3 ^/D 6 
Using this matrix one can write all polarization observables in the same form as for 
the elastic channel. In particular the final proton analasing power will be given by the 
formula 

Trqt^u } 

Ti{MA} V ; 

where 0V M1 is the spin matrix of the 1st nucleon, supposed to be proton (n = +1/2). 

3. Results 

We present our results for the differential cross section, nucleon analyzing power A y , 
and deuteron vector iT n for nd elastic scattering at 14.1 MeV in Figs. [TH3]In all our 
calculations we used the AV14 NN potential. The total angular momentum of the pair 
of nucleons j'23 has been taken up to 3 and all values of the total three-nucleon angular 
momentum M up to 13/2 with both parity values (±1) have been taken into account 
(up to 62 partial waves in all). 

Fig. [1] shows our results for elastic differential cross sections along with the 
prediction of the Bochum group Ref. [21] using the charge-dependent AV18 NN 
potential. In Fig. |2] our results for nucleon analyzing power are shown together with 
the prediction of the Bochum group and experimental data from [22]. In general 
the agreement between our calculations of A y and those of the Bochum group is 
satisfactory, since using different NN potentials may explain small differences between 
these calculations around the maximum values of A y (Fig. [2]) 

Our predictions for the deuteron analyzing power iT n and those of the Grenoble 
group [23] are presented in Fig. [3] along with the experimental data [24J . Again there is 
a small disagreement between two calculations which can be explained by the use of two 
different type of NN potentials. The results of the Grenoble group have been obtained 
with NN interactions of de Tourreil and Sprung, effective in the states ^o, 3 S\, 3 Di, 1 P±, 
and 3 P ,i,2- 

For nd breakup scattering at 14.1 MeV, our preliminary results for the breakup 
differential cross section and nucleon analyzing power A y under FSI configuration are 
shown in Figs. I4"f31 along with the experimental data [23] . Results for A y are presented 
together with the calculations of Witala et al. [26] using the Paris NN potential and 
Bonn potential of version A and B with total two-nucleon angular momentum j '< 2. 
Thus in their work up to 34 partial wave states have been taken into account for each 
total angular momentum J and parity states. 
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4. Discussion 

Our results for nd elastic scattering at 14.1 MeV and those from the Bochum and 
Grenoble group are in fair agreement. Some differences can be attributed to smaller 
values of the total two-nucleon angular momentum j taken into account in their 
calculation as well as to different NN potentials used. In the energy region of tens MeV 
[22] and [24] theoretical predictions are 25-30% lower than the experimental data. 

For nd breakup observables our predictions compare only qualitatively with the 
experimental data and predictions for A y of the Bochum group. Still the positions of the 
peaks in the breakup angular distribution are correct. For the nucleon analyzing power 
the deep minimum at arc length S equal to 7 in our results disagrees with predictions 
for A y from [26J. 

5. Conclusion 

For nd elastic observables at 14.1 MeV acceptable agreement (in view of different 
NN potentials) between predictions of our calculations and those of the Bochum and 
Grenoble group demonstrates the soundness of our novel method providing thereby a 
new approach for calculating three-nucleon scattering. Our approach can and will be 
used to include the Coulomb force. It is well-known that Nd polarization observables 
are the magnifying glass for studying 3 Pj NN states and calculations which rigorously 
include nuclear and electromagnetic interactions are very important. 

Our next step is to study pd breakup scattering using a new and more accurate 
version of the Faddeev equations in configuration space. In presence of the Coulomb 
interaction serious changes are to be made in view of p and n being different particles, 
so that use of the isotopic formalism becomes unreasonable. As a result the Faddeev 
equations in configuration space derived by S.P. Merkuriev et al. in [17] become 
only approximate. Probably the error of their calculations of the elastic pd scattering 
amplitudes is not significant. However in the case of pd breakup scattering one has to 
use both charge-dependent NN AV18 potential and the new correct Faddeev equations. 
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Figure 1. The differential cross section for nd scattering at Ei a b=14.1 MeV. The solid 
line is our results. The dashed line corresponds to prediction of the Bochum group 
Ref. ED. 
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Figure 2. The nucleon analyzing power A y for nd scattering at Ei a b =14.1 MeV. The 
solid line is our results. The dashed one is results of the Bochum group Ref. [21.. The 
experimental data are from Ref. [22] . 
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Figure 3. The deuteron analyzing power iTn for nd scattering at Ei a b =14.1 MeV. 
The dot-dashed line is prediction of the Grenoble Group with SSC NN interaction Ref. 
[23] . The proton-deuteron data at Ei a b =15.0 MeV are from Ref. [24] . 
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Figure 4. nd breakup differential cross section as a function of the arc length S 
under FSI configuration (6>i =52.6°, 6 2 = 40.5°, <j> 12 =180°) at E lab =14.1 MeV. The 
proton-deuteron experimental data denoted by solid triangles are from Ref. 25 . 
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Figure 5. The nucleon analyzing power as a function of the arc length S under FSI 
configuration (6 X =52.6°, 6 2 = 40.5°, <f> 12 =180°) at E lab =14.1 MeV. The solid line is 
our results. The dashed one is results of Ref. [26 . The proton-deuteron experimental 
data denoted by open squares are from Ref. [25]. 



